Towards an experimentally feasible controlled-phase gate on two blockaded Rydberg 

atoms 



Michael Murphy0 Simone Montangero, and Tommaso Calarco 
Institut fiir Quanteninformationsverarbeitung, Universitdt Ulm, 
Albert- Einstein- Allee 11, D-89069 Ulm, Germany 



o 

(N 

> 
O 

in 



> 
cn 

00 

o 



Philippe Grangier and Antoine Browaeys 
Lahomtoire Charles Fabry, Institut d'Optique, CNRS, Univ Paris-Sud, 
Campus Poly technique, RD 128, 91127 Palaiseau cedex, France 

We investigate the implementation of a controUed-Z gate on a pair of Rydberg atoms in spatially 
separated dipole traps where the joint excitation of both atoms into the Rydberg level is strongly 
suppressed (the Rydberg blockade). We follow the adiabatic gate scheme of Jaksch et al. where 
the pair of atoms are coherently excited using lasers, and apply it to the experimental setup outlined 
in Gaetan et al. @]. We apply optimisation to the experimental parameters to improve gate fidelity, 
and consider the impact of several experimental constraints on the gate success. 



I. INTRODUCTION 

Using neutral atoms for quantum information has gar- 
nered much theoretical interest over the last decade, fu- 
elled by advances in their experimental manipulation, 
particularly trapping and cooling. Several novel schemes 
for entangling pairs of atoms (an essential operation for 
quantum logic) via controlled collisions have been devel- 
oped I, g, but schemes that make use of the special 
properties of Rydberg atoms are also very promising (see 
[5| for a review). In particular, several schemes for pro- 
ducing quantum gates by excitingpairs of Rydberg atoms 
with tuned lasers have emerged |6| which capitalise on 
the strong dipole-dipolc interaction that prevents the si- 
multaneous excitation of neighbouring Rydberg atoms, 
known as the Rydberg blockade. Several steps towards 
realising such schemes experimentally have already been 
achieved, particularly the observation of the blockade [ll 
and entanglement generation Q in a system of two con- 
fined Rydberg atoms. There has even been some early 
success in producing a gate with trapped Rydberg atoms 



the movement of the atoms in the dipole trap. We will 
attempt to address both of these issues here by apply- 
ing a gradient ascent control algorithm to search for the 
ideal set of parameters for implementing the gate on a 
short timescale 1 fis) and with high fidelity. Our re- 
sults will show that the physical system allows for a great 
deal of control and gate times and fidelities approaching 
our desired range, providing a positive outlook for imple- 
menting high-fidelity gates with such systems. 

The paper is arranged as follows. In Sec. |TT1 we briefly 
recount the CZ gate, followed by a description of how 
it may be synthesised on a pair of Rydberg atoms. In 
particular, we expose the operation of the gate by con- 
sidering the effective two-level dynamics of each atom 
and its interaction with a laser field. In Sec. IIIH we de- 
scribe how we optimise the operation of the laser using 
a gradient descent to achieve the gate with high-fidelity. 
In Sec. IIVI we consider the details of the experiment and 
the constraints it imposes on the gate operation, particu- 
larly with regards to loss and effects arising from atomic 
motion. Finally, we conclude our paper in Sec. |Vl 
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In this paper, we consider the implementation of a 
controlled-Z (cz) gate on a pair of Rydberg atoms con- 
fined in spatially separated dipole traps subject to the 
Rydberg blockade effect. We follow the scheme outlined 
in but with specific application to the experimental 
setup detailed in Q, where the Rydberg atom is excited 
via a two-photon transition. This proposal has the ad- 
vantage that both atoms are excited by the same laser, 
reducing the need for single-atom addressability; the gate 
is also adiabatic, which softens the experimental require- 
ment for strong fields or precise timings. However, the 
experimental considerations do present additional chal- 
lenges in the implementation of the gate, particularly due 
to loss from the intermediate state of the transition and 
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II. CONTROLLED PHASE GATE 

A. Gate definition 

The controlled-Z (cz) gate is a two-qubit gate in quan- 
tum information, and belongs to the class of controlled 
unitary operations [loj . Given the computational basis 
|0), |1), it is defined as the unitary transformation 
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This gate is of particular importance because it generates 
entanglement between two unentangied qubits. In addi- 
tion, together with a finite set of single-qubit operations, 
one can construct any desired gate operation simply by 



2 



S{ 



1.9) 



FIG. 1. The level scheme for a single Rydberg atom. The 
atom is driven to the Rydberg state via a two-photon transi- 
tion, which couples the ground state \g) to the excited state 
\r) through the intermediate state The effective Rabi fre- 
quencies of the transitions are finXi) for the red laser (which 
is blue-detuned by A) and Qb for the blue laser (which is 
red-detuned by A -f 5). 



taking combinations of these operations with the CZ gate. 
This is known as a universal set for quantum computa- 
tion [Til. 



in several papers [J H, [l^ [l^l, but here we follow the 
adiabatic (model B) scheme of Jaksch et ai, where the 
gate is performed by adiabatically driving the two-atom 
system 

C. Gate operation in outline 

There are two critical elements that allow us to sythe- 
sise the gate with our system. The first is the block- 
ade mechanism, which prevents excitation to the doubly- 
excited \rr) state (where we have used the shorthand no- 
tation I r) ® I r) — I rr) for the tensor product of the state of 
the two atoms, which will used throughout). This avoids 
unwanted mechanical effects stemming from the strong 
interaction of the two Rydberg atoms, as well as reduc- 
ing the time spent in the Rydberg state, which is subject 
to loss. The second crucial aspect is the super-radiant 
enhancement of excitation from \gg) as compared with 
the states \ge) and \eg). This results in a higher rate of 
phase accumulation on the \gg) state during excitation 
in comparison to \ge) and \eg). By carefully choosing the 
excitation profile of the incident lasers, we can control 
these two different accumulated phases to produce the 
CZ gate. 



B. Blockaded Rydberg atoms 

The physical system we are considering for the imple- 
mentation of the gate is a pair of trapped Rydberg atoms 
0, specifically ^^Rb. The atoms are trapped a distance 
r apart in two separate microscopic dipole traps Q . For 
our purposes, we need only consider a small number of 
the internal states on which the dynamics will take place. 

A pair of hyperfine states of the atom will encode the 
computational basis, and are labelled |e) = |0), \g) = 
Each atom has its \g) state coupled to a highly excited 
Rydberg state (which we label |r)) via a two- photon tran- 
sition through an intermediate state \i). The internal 
level scheme with state transitions for a single atom is 
shown in Fig. [TJ If two neighbouring atoms are excited to 
the \r) state, then they interact via a dipole-dipole poten- 
tial with energy U{r) = C^/r^. This interaction energy 
shifts the energy of the state where both atoms are ex- 
cited. When this shift is much larger than the two-photon 
detuning S, the two-photon transition is far off-resonant 
with the doubly-excited state, leading to a strong sup- 
pression of both atoms becoming excited. This effect is 
known as the Rydberg blockade, and has been observed 
experimentally [2, [13 ■ The effect is shown schematically 
in Fig. [2 

The important point about the blockade mechanism 
in our case is that it is state dependent: only if both 
atoms are in the ground state \g) will they be subject 
to the blockade. The potential use of this as a mecha- 
nism for performing a quantum gate has been explored 



D. Hamiltonian 

The two-photon transition is driven via two lasers; 
one blue-detuned on the transition from \g) to \i) by an 
amount A with a Rabi frequency flji(t), and the other 
red-detuned on the transition |i) to |r) by an amount 
A + S with a Rabi frequency Hb (sec Fig. [T]). The gen- 
eral form of the Hamiltonian for our two-atom system 
can be written as 



H — H,. + + Hi, 



(2) 



The single-atom Hamiltonians are composed of both the 
internal and external dynamics, such that (after the ro- 
tating wave approximation) 



{A-i-f,)\i){i\ + {d -i-fr)\r){r\ 

fif^fl(0(^.(k.r,+At)|g)^,H_H.c)- 
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i(k-r3+(A+<5)t) 



(3) 
(4) 

(5) 



i)(r| +H.c) 



Hi = (f + ytrap)(|5)(.g| + |e)(e| + \^){^\ + \r){r\) , (6) 

where i — 1,2 labels the two atoms. Hs describes the 
energy splitting of the internal states along with the ef- 
fective decay from those states, Hj describes the laser 
coupling between the internal states, and He contains 
the kinetic and potential energy terms. The factors 7^, 
7r account for an effective loss of population from the 
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FIG. 2. The level scheme for the two-atom system (neglecting 
the internal level The transition from the joint ground 
state \gg) to the super-radiant state — {\gr) + \rg))/^/2 

is enhanced by a factor of \/2, while there is no coupling to 
the sub-radiant l^*") = (Igr) — |rgi))/\/2 state. The excitation 
of both atoms to the Rydberg state \rr) is forbidden, since the 
interaction energy U{r) has shifted the level far off-resonant 
with the incident lasers. 



intermediate and Rydberg levels respectively. The ex- 
ponential terms in Eq. ([S|) arc phases accumulated by 
the atoms as they move through the light-field of the 
laser; k is the wavcvcctor of the laser field, and rj is the 
position vector of the jth atom. The interaction Hamil- 
tonian is given simply by the dipole-dipole interaction: 
Hint = U{r)\rr){rr\. 

Note that in Eq. (jB]) we have neglected the difference 
in the trapping potentials for the different internal states 
(in any case we will turn off the trap when the gate is 
performed). Since there are no state-dependent terms in 
Eq. ([S]), we can neglect it in our treatment. 



We can now rewrite the system Hamiltonian as H 
Hr + Hint , where 



Hr = + + Hi 



Hi 



^^''%){9\ + {5-'^-^lr)\r){r\ 
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(8) 
(9) 



Hi = -^(e-"'--Mge)(re| +e-*--^|eg)(er| +H.c.)- 



V2^(|.g.g)(*+|+|*+)(.g.g|) 



(10) 



with the effective Rabi frequency of the two-level dynam- 
ics 

n{t) ^nB^Rit)/2A. (11) 

The magnitude of the dipole matrix elements 



\{gg\Hi\i^+)\ ^ V2 ■ 



\{gg\Hi\^p-)\=0.^ 



\{ge\Hi\re)\ - \{eg\Hi\er) 



show that the state l^*^) is not coupled to any of the 
other states via the laser interaction; in other words, it is 
sub-radiant In addition, the state l^""^) is a super- 
radiant state, so that the coupling between the ground 
state \gg) and 1^'"*') is enhanced by a factor of -\/2 com- 
pared to the transition \ge) — >■ |re) {\eg) — > |er)). 

Finally, note that since \i) is never populated in this 
approximation, so we neglect the loss term 7^. We will 
also for the moment assume that the atoms are station- 
ary, and so the phases accumulated from their movement 
in the light-field can be neglected. 



F. Gate operation in full 



E. Effective two-level system dynamics 

By performing an adiabatic elimination [l5l | of the 
state we can examine the effective three-level dy- 
namics of the system. This leads to the condition that 
{lR{t),i}R,S <C A. In addition, we make a change of 
basis such that 

\gr) ^ |vl/+) ^ ± (e^k.r,|^^^ ^g.k.r.|^^^) ^ 

\rg) ^ I*-) ^ ^ (e^-^-"-^ \gr) - e^^^^^ \rg)) . (7) 



Now we describe the operation of the gate in more 
detail. We start with the initial state 

Hit = 0)) = i i\gg) + \ge} + \eg) + \ee)) . (12) 
and define the target state at final time T 

Hg) = I Hgg) + \ge) + \eg) + \ee)) . (13) 

Note that while this seems to a specific state transfor- 
mation, as opposed to the unitary transformation from 
Eq. ([T]), they are in this case equivalent by virtue of the 
basis states {\gg), \ge)^ l^ff), |ee)} not being directly cou- 
pled to one another. Hence any initial state-dependent 
phases will not affect the final outcome of the gate. 

To perform the gate, the blue laser is always switched 
on, while the red laser is modulated in a time-dependent 
fashion using an acoustic-optical modulator. If this mod- 
ulation is slow on the timcscale given by n{t) and 5, then 
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the system will adiabatically follow the dressed states of 
the Hamiltonian H. Performing the same treatment as 
in [l| for our system, we similarly find that the energy 
of the dressed levels adiabatically connected to \gg) and 
l5e)(|e5)) are 



-ge 



6' - 2ERit) - 



(14) 
(15) 



respectively, where 6' = 5 — Eb + Eji{t) is the effective 
two-photon detuning including the Stark shifts from the 
adiabatic elimination of \i): 



ER{t) 



4A 



and 



Ef 



4A 



2u + 46' - iEnit) 



(16) 



(17) 



includes the additional Stark shift from the adiabatic 
elimination of the \rr) state. The final state is 



where 



'Igg) 



{\ge) + \eg)) + \ee), (18) 



'99 



(t) dt , 0, 



9e 



-ge 



it) dt , (19) 



By performing state-selective qubit operations, we can 
realise the cz gate. To see this, we first apply a state- 
selective phase on the first atom: if the first atom is in 
the state \g), then it receives a phase 6""^'= . Similarly, we 
then apply the same phase rotation on the second atom. 
After these operations, the state becomes 



|^(r)) = e-'(*--2*-)|5g) + |.ge) + leg) 
We can now define the gate phase 
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(20) 



(21) 



The operation of the gate is now clear: we seek to modu- 
late ftR{t) such that (f) ~ (2fc + l)7r, k £ Z, and there is no 
remaining population in the excited states of either atom 
(this is taken for granted in the adiabatic limit). The 
next step is to design ilR{t) to achieve these conditions. 



III. OPTIMISATION OF THE GATE 

A. Simulation 

The system evolves in accordance with the Schrodinger 
equation (h = 1): 



(22) 



Initial 300 300 1000 

Optimised 304.66 292.55 974.78 

TABLE I. A set of initial parameters that produces a gate 
with fidelity of around 89%, and the optimised parameters 
that produce a gate with fidelity better than 99.9%. All values 
are in units of 2n MHz 



Since the Hilbert space dimension \H\ = N is relatively 
small, we can simulate the gate by directly diagonalising 
the Hamiltonian and using discrete time steps dt, such 
that the solution of Eq. ([22]) can be written as 



\m)^pe-''"''p-'mt-dt)) , 



(23) 



where P — [xi,X2t ■ ■ iXn] is the square matrix con- 
structed from the eigenvectors x,; of H, and D = 
diag(Ai, A2, . . . , Aat) is a diagonal matrix whose elements 
are the eigenvalues of H. 



B. Optimisation method 

Many tools exist for numerical optimisation that we 
could employ here to design the Rabi frequency fi_R(t) 
that produces the gate. We start by making a guess for 
the form of ilR{t), so that the evolution is characterised 
by only a handful of parameters. One particular choice 
is a Gaussian: 



Ofl(t) = n 



e 



(24) 



where we will choose r = 0.1. The constant Hq is a pa- 
rameter which we can, in principle, choose arbitrarily (in 
reality there will be constraints on this value, which we 
will come to later) . There are also parameters associated 
with the system that one may vary, namely the atom 
separation r. the Rabi frequency fls, and the detunings 
S and A. We begin by choosing a set of reasonable val- 
ues, and then we numerically optimise each of the pa- 
rameters to achieve the desired gate with a high fidelity. 
The numerical method used to optimise the parameters 
is gradient descent. 

As an example, we start with the set of parameters in 
TableU The total time for the gate is fixed at 500 ns, and 
the atoms are separated at a distance of 0.3 /.tm, and the 
detuning 6 = 27r300 MHz. The dipole-dipole interaction 
strength U{r) = 3200MHz/r3 is taken from 0. This is 
a particularly strong interaction due to the existence of 
a Forster resonance [itI . [18| , a point which we will come 
back to later when we describe the experimental setup in 
more detail. We define the fidelity of the gate operation 
as F = |(V'(T)|V'g)P- The resuhant fidelity of the gate 
with these initial parameters is only around 89%. Now we 
apply our optimisation algorithm to the parameters Vlo, 
Hb, and A. After 1600 iterations, we achieve a fidelity 
of better than 99.9%, or, more precisely, and infidelity 
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FIG. 3. The decrease in infidelity of the quantum gate for the 
set of initial parameters given in Table H] One sees that the 
infidelity decreases monotonically until saturating in a local 
minimum of the optimisation. The final achieved infidelity 
was 1 - = 3.8 • 10"'*. 
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FIG. 4. (Color online) The dashed (green) line is the phase 
accumulation of the states \gg) (given by (pgg), the dotted 
(blue) line that of \ge) (given by (jige), and the solid (red) line 
is the total entanglement phase 0, which at the final time 
reaches 3.0 (note that the dashed green line and the solid red 
line overlap almost exactly). The parameters used are given 
m Tab. m The time is in units of ns, and so the duration of 
the gate is 500 ns. 



/ = (1 — F) < 4 • 10"'^. The decrease in infidelity as the 
algorithm progresses is shown in Fig. |31 Table |T] shows 
the final set of parameters that produced the optimised 
gate. The resulting time-dependent phase accumulations 



and 4>ge are shown in Fig. |31 We see that the phase 
accumulated by the £ge{£eg) is exactly zero. This can be 
understood by calculating Eq. ([T5|) for our set of param- 
eters: since i7o = ^^Si it is straightforward to show that 
£ge{t) = 0. The population of the relevant levels during 
the gate operation are shown in Fig. [5] 

It might be worrisome that the parameters in Tab.|T]do 
not seem to strongly fulfil the condition f2ij(i), fiij, (5 ^ 
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FIG. 5. (Color online) (a) The solid (red) line is the pop- 
ulation of the state \gg) and the dashed (green) line is the 
population of the state 1^'''") over the duration of the gate, 
(b) Similarly to (a), the solid (red) line is the population of 
the state \ge) and the dashed (green) line is the population of 
the state \gr) over the duration of the gate. The populations 
of \eg) and \er) are respectively the same. 



A required for the validity of the adiabatic elimination. 
We have, however, confirmed that even when considering 
the full evolution under H in Eq. ^ the population in |i) 
is heavily suppressed, such that it may be neglected. In 
what follows, we will abandon the effective model given 
by Hji (which provided insight into the gate mechanism) 
in favour of the full treatment by H. 



IV. EXPERIMENTAL CONSIDERATIONS 



A. Level description 

We have now demonstrated the optimisation method, 
but we must also consider the experimental conditions 
which will have an effect on the gate fidelity. We con- 
sider the setup given in The gate is well suited to 
this system because we do not require single atoms to 
be addressable, and the Rydberg blockade is relatively 
strong. The reason for this is the use of a Forster res- 
onance that exists in ®^Rb which comes about due 
to the quasi-degeneracy of the two-atom states (58c?3/2, 
58^3/2) and (60pi/2, 56/5/2)- This enhances the dipole 
interaction, leading to an interaction energy U{r) oc Ijr^ . 

The choice of states for the levels are |g) = |5si/2, F = 
1, Mf = 1), |e) = |5si/2, F = 2, Mp = 2), and \i) = 
|5pi/2, F ~ 2, Mp — 2), as in reference 0. 



B. Correspondence to experimental results 

As mentioned earlier, we describe the loss from both 
the intermediate state \i) and the excited Rydberg state 
I r) phenomenologically through the decay rates ji and 7^ 
respectively. Examining the literature [3, [l^ , we find for 
our setup that ji ~ 2tt ■ 5.75 MHz and 7^ = 27r • 4.8 kHz. 
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FIG. 6. (Color online) In (a), the solid (red) line is the simu- 
lated probability to excite \ge) to |re), with the experimental 
results for the same transition shown in dashed (green). Sim- 
ilarly in (b), the solid (red) line is the simulated probability 
to excite \gg) to |^^), with the experimental results for the 
same transition shown in dashed (green). 



We have verified that the simulation corresponds to 
the experimental results from Q. Fig. [S] shows the 
agreement between the simulated Rabi oscillations of the 
\ (|er) + |re)) state and the state. While the fit 

is not exact, we do reproduce the correct frequency of 
oscillation, as well as an indication of the typical decay 
from the intermediate level. We also find a relative dif- 
ference in frequency of the two oscillations of a factor 
~ \/2, as expected from the theory. The discrepancies 
between theory and experiment arise from experimental 
imperfections which are not taken into account in our 
model, namely laser fluctuations in both power and fre- 
quency of the lasers, which lead to some dephasing. Our 
phenomenological loss model also does not account for 
the possibility that an atom can decay to \g) from where 
it may be repumped, which partially accounts for the 
discrepancy in total population. 



expect that the final result will not depend very strongly 
on the exact shape of the excitation as long as the area of 
the pulse is preserved.) The 'rise-time' (the time it takes 
to increase the laser power from zero to its maximum) 
is typically in the range 200-400 ns. We take the pro- 
file of this rise-time to be Gaussian. Lastly, due to the 
separate dipolc traps, the minimum distance between the 
two Rydbcrg atoms is limited to 3 jiva or above. A fi- 
nal condition is put on the detuning A of the red laser: 
this should be less than 500 MHz due to experimental 
constraints (although this is not a stringent condition). 

There is also an additional probability of loss when the 
atoms is excited to the Rydberg level: when excited, the 
motional wavepacket starts to spread, so that when the 
dipolc trap is reapplied, there is a finite probability that 
the atom is lost (this is actually used as a method of 
detection in the experiment). This motivates us to limit 
the time spent in the Rydberg state. 



D. Optimising the gate for experiment 

With these limitations, we now see that the gate pa- 
rameters from Tab. [T] are not feasible in our chosen ex- 
perimental setup. We must now start with a new set 
of parameters and run the optimisation again to see to 
what extent the gate is still implcmcntable. Given the 
discussion above, we are motivated to make the following 
changes to our parameters. 

• The gate should be performed as quickly as pos- 
sible, meaning that the effective Rabi frequency 
VL (X Ofl, JIb should be made large. This implies 
that we should choose 51/? and fls close to their 
maximum values. (This has the additional advan- 
tage that we spend less time in excited Rydberg 
levels, improving the probability of recapture.) 

• To avoid excitation of the lossy \i) state, we need 
to keep the red laser far-detuned, ideally around 
500 MHz. However, as can be seen in Eq. (|lip . 
increasing the detuning will reduce the effective 
Rabi frequency, which makes achieving the gate in 
a short time more difficult. 



C. Typical experimental parameters 

While it would be ideal if the gate parameters de- 
scribed in the last section could be immediately applied 
in the experiment, the reality is that there arc certain 
experimental limitations that prevent us from doing so. 
Firstly, the power of both lasers have certain maximum 
values: flu has a maximum operating value of 400 MHz, 
while Hb is limited to 30 MHz. Secondly, we do not 
have the freedom to modulate the laser power as we like; 
acoustic-optical modulators (AOMs) control the power 
of the laser beams incident on the atoms, and they have 
limits on the rate at which the intensity of the beam can 
be changed. (In any case, since the gate is adiabatic, we 



• To make maximum use of the laser power, we 
change the shape of the pulse from a simple Gaus- 
sian to a 'flat-top', given by 




t+sgn(f)(r-T/2)^ 
7^78 



|i|>i-' 

otherwise. 



(25) 



Here, the Rabi frequency increases with a Gaus- 
sian profile to the maximum in a time r (the 
rise-time), followed by a period of constant Rabi 
frequency for a time T — 2r, and then finally a re- 
duction along a Gaussian profile to zero, again in 
a time r. This allows us to have the laser at full 
power for the longest time possible, which in turn 
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flo A 5 

Initial 49 28 500 

Optimised 52.31 26.85 496.81 -1.12 

TABLE II. A set of initial parameters, within experimental 
constraints, that produces a gate with fidelity of around 32%, 
and the optimised parameters that produce a gate with fi- 
delity better than 96%. All values are in units of 2n MHz. 

1 P < 1 




10-2 I , 1 

500 1000 

Iterations 

FIG. 7. The decrease in infidelity of the quantum gate for the 
set of initial parameters given in Table [III One sees that the 
infidelity decreases monotonically until saturating in a local 
minimum of the optimisation. The final achieved infidelity 
was 1 - F = 3.1 • 10"^ 

causes us to accumulate the time-dependent phase 
more rapidly. 

• Since the minimum distance between the atoms in 
the experiment is r^in — 3.0 /im, we will use this 
value in what follows to ensure we are as deep in 
the blockade regime as possible. 

Based on these considerations, we try the set of pa- 
rameters given in Tab.HIl As before, we can examine the 
convergence of the optimisation (Fig. [7]), the accumula- 
tion of the entanglement phase (Fig. [5]), and the shape 
of the Rabi frequency ilR{t) (Fig. |9]). The final gate fi- 
delity is only around 96%, but we are able to successfully 
avoid losses, resulting in a final population norm of 0.995. 
This high avoidance of loss can be seen in the high con- 
trast Rabi oscillations in Fig. [TOl By comparison with 
the results in Fig. [TTJ we can also see a relative differ- 
ence in Rabi frequency of ^ between the transitions 
\gg} 1 4'+) and \eg) ^ |er). 

We believe that this result is close to the optimal case 
for this system and gate implementation given the exper- 
imental limitations. There are essentially two competing 
factors which makes performing the gate more difficult; 
one the one hand, the effective Rabi frequency should 
be made large to increase the rate at which we accu- 
mulate the phase, but this is countered by the need to 
make the single-photo detuning A as large as possible to 
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FIG. 8. (Color online) The dashed (green) line is the phase 
accumulation of the states \gg} (given by (pgg), the dotted 
(blue) line that of \ge) (given by 4>ge), and the solid (red) line 
is the total entanglement phase (j>, which at the final time 
reaches 1.0. The parameters used are given in Tab. [Ill The 
time is in units of ns, and so the duration of the gate is 2 fis. 
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FIG. 9. (Color online) The pulse shape for the Rabi frequency 
n_R(t) for the parameters in Tab.HTl The shape is a 'fiat-top', 
as given in Eq. (f25|l . 



avoid populating the intermediate state \i). In addition, 
this inevitably leads to long gate times on the order of 
microseconds, which increases the chances of losing the 
atom while in the Rydberg state (we approximate that 
each atom spends 500 ns in the Rydberg state for the 
gate with duration 2 /us). 



E. Movement of the atoms in the hght-field 

Until now, we assumed that the atoms were stationary 
in the dipole traps. In reality, the atoms are Doppler- 
cooled to around 75 fiK. During the laser excitation, the 
trapping fields are switched off, allowing the atoms to 
move freely in any direction in the plane perpendicular 



8 




I > > 1 1 
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t (ns) 

FIG. 10. (Color online) The solid (red) line is the population 
of the state \gg) and the dashed (green) line is the population 
of the super-radiant state over the duration of the gate. 




-1000 -500 500 1000 

t (ns) 

FIG. 11. (Color online) The solid (red) line is the population 
of the state \ge} and the dashed (green) line is the population 
of the state \gr) over the duration of the gate. The popula- 
tions of \eg) and \er) are respectively the same. 

to the trapping field. The terms -di = arccos(k • r^) e 
[0,7r] from Eq. ([TU| then produce additional independent 
phases on each atom. Since we don't know a priori in 
which direction the atoms will move with respect the 
light field, the phase difference between the two atoms is 
essentially random. 

Figures [1^ and [T2] show the effect of this phase on the 
final fidelity of the gate for the parameters for the sets of 
optimal parameters in Tabs. HI and HIl respectively. In the 
first implementation of the gate, we see that the effect 
of the motional phase on the gate fidelity is, at worst, 
a drop in infidelity from 10"'' to 10~^. In the second 
case, the effect is much more significant, in some cases 
producing fidelities near zero. This is mostly due to the 
increase in the gate time T. 

To actually perform the gate in practice under these 
conditions, we would have to cool the atoms much closer 




FIG. 12. (Color online) The effect of the motional phases 
of the atoms on the final gate fidelity for the parameters in 
Tab. in Here, the effect of the motional phases is not too large, 
since we perform the gate over a very short time and we do 
not significantly excite the atoms. 




n/2 



FIG. 13. (Color online) The effect of the motional phases 
of the atoms on the final gate fidelity for the parameters in 
Tab. [nl Here, the effect of the motional phases is detrimental 
to the gate fidelity. 

to the motional ground state. This would reduce the 
distance that the atoms move in the light-field, and hence 
the amount of phase that they collect. The dependence 
of the fidelity on the average temperature of the atoms is 
given in Fig 1141 Thus by reducing the temperature by a 
factor of around five will increase the fidelity to around 
90%. This is experimentally realistic as demonstrated in 
reference [1^. 

V. CONCLUSION 

We have investigated the implementation of the Ryd- 
berg two-qubit entangling gate from Jaksch et a/.[l| in 
the experiment outlined in We applied a gradient 
ascent control algorithm to find optimal sets of exper- 
imental parameters that produced a gate with low loss 
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Temperature {fJ.K) 

FIG. 14. (Color online) The final fidelity of the gate plotted 
against the average temperature of the atoms. Each (green) 
point is an average over 10,000 realisations of i9i and chosen 
randomly in the range [0,7r]. The (red) solid line is a linear 
fit to the points. The (blue) dashed line shows the current 
temperature of 75 /iK. 



from the intermediate level while still achieving a fidelity 
of around 96%. The main source of error was found to 



be the random phase accumulated by the atoms in the 
light field. 

While the system seems ideally suited for this scheme, 
the experimental limitations do not allow for a great deal 
of success in implementing the gate. The most notable 
source of error comes from the movement of the atoms 
in the light field, which could be minimised in future 
experiments by cooling the atoms further. Eliminating 
this source of error should allow gate fidelities of around 
99% which, while not quite good enough to allow quan- 
tum computation (even with error correction), would be 
a significant step forward for the realisation of quantum 
computation with neutral atom systems. 

It is worth pointing out that while we investigated the 
gate scheme of there are some alternative schemes 
that could be implemented in our setup, most notably 
perhaps the scheme of @ which uses a STIRAP pulse se- 
quence to excite the atoms. We have also not investi- 
gated allowing the Rabi frequency of the blue laser fis 
to modulate in time which could lead to a more robust 
gate implementation. 

We acknowledge financial support by the EU un- 
der the contracts MRTN-CT-2006-035369 (EMALI), IP- 
AQUTE, and from the German SFB TRR21. We also ac- 
knowledge computational resource provided by the BW- 
grid. 
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